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Executive  Summary 

The  National  Search  and  Rescue  (SAR)  Manual  describes  a  method  for 
calculating  wind  driven  currents  which  is  based  on  a  manual  method  devel¬ 
oped  by  Mooney  (1978).  This  report  discusses  a  preliminary  investigation  of 
the  model  to  determine  its  usefulness  in  search  and  rescue  applications.  The 
derivation  of  the  model  has  been  reviewed  to  verify  that  it  is  physically  cor¬ 
rect.  It  has  been  compared  with  the  Mellor-Yamada  Level  2|  model  (Mellor 
and  Yamada,  1982),  a  complex  one  dimensional  model  of  the  mixed  layer. 
Several  aspects  of  the  model  have  been  assessed  and  the  conclusions  will  be 
summarized  below. 

In  comparison  with  the  Mellor-Yamada  Level  model,  the  SAR-Mooney 
model  over-estimated  the  deviation  of  the  current  direction  from  the  wind 
direction  for  the  entire  range  of  wind  velocities  examined  (2.5  to  40  m/s  [5 
to  80  knots]).  The  model  predictions  of  current  velocity  agreed  well  with 
a  steady  wind  speed  of  10  m/s  (20  knots),  but  the  difference  between  the 
models  increased  with  increasing  windspeed. 

The  SAR-Mooney  model  was  tuned  with  subsurface  current  data  from 
12  m,  however,  when  compared  with  the  Mellor-Yamada  model,  it  agreed 
with  near-surface  predictions.  Although  this  result  suggests  that  the  SAR- 
Mooney  model  may  not  describe  well  the  experimental  data  on  which  its 
coefficients  are  based,  it  does  not  appear  to  be  a  drawback  to  its  use. 

Mooney  (1978)  determined  his  model’s  coefficients  with  data  from  rela¬ 
tively  deep  water  compared  with  the  depth  limitation  recommended  by  the 
National  SAR  Manual  of  30  m  (100  feet).  Where  the  water  depth  is  equal 
to  or  greater  than  the  Ekman  depth  which  is  the  depth  of  the  bottom  layer 
influenced  by  the  wind  stress,  the  wind  drift  should  not  be  affected  by  the 
water  depth.  At  depths  shallower  than  the  Ekman  depth,  the  SAR-Mooney 
predictions  should  be  treated  with  caution.  Figure  1  may  be  used  to  estimate 
the  limiting  depth  for  which  the  water  can  be  assumed  to  be  infinitely  deep 
as  a  function  of  windspeed  and  latitude. 

Mooney  (1978)  tuned  his  model  with  wind  and  current  data  from  about 
39°  N.  It  appears  that  the  SAR-Mooney  model  may  be  increasingly  inaccu¬ 
rate  as  the  equator  is  approached,  particularly  with  respect  to  the  magnitude 
of  the  current. 

In  the  vicinity  of  the  equator,  the  inertial  frequency  approaches  zero,  so 
the  wind  drift  is  dependent  only  on  the  wind  stress  and  therefore,  the  drift 
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will  be  in  the  direction  of  the  wind. 

The  recommended  limitation  of  the  SAR-Mooney  method  to  cases  greater 
than  32  km  (20  miles)  from  the  coast  appears  to  be  too  short  a  distance. 
Wind  drift  will  be  affected  by  a  coastal  barrier  at  a  distance  equal  to  or 
less  than  the  Rossby  radius.  As  the  water  depth  decreases,  the  effect  of  the 
coastline  increases.  Additionally,  the  effect  of  the  coast  increases  linearly 
with  decreasing  distance  from  the  coast. 
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1  Introduction 


A  quantitative  understanding  of  the  mechanics  of  wind  generated  currents 
is  of  great  importance  to  oceanographers  and  others  concerned  with  such 
practical  problems  as  search  and  rescue,  oil  spill  propagation  and  iceberg 
drift.  These  types  of  problems  have  a  need  for  near  real-time  estimates 
of  the  local  wind  driven  current.  The  National  Search  and  Rescue  (SAR) 
Manual  describes  a  technique  for  calculating  the  local  wind  current  which  is 
based  on  a  manual  method  developed  by  Mooney  (1978).  This  method  has 
also  been  incorporated  into  Computer  Assisted  Search  Planning  (CASP)  and 
the  iceberg  drift  model  used  by  the  International  Ice  Patrol  (IIP). 

The  purpose  of  this  report  is  to  examine  the  National  SAR  Manual 
method  for  determining  local  wind  driven  current  and  to  address  the  fol¬ 
lowing  questions:  (1)  is  the  SAR  method  physically  correct;  (2)  does  the 
SAR  method  accurately  describe  the  wind  generated  currents;  (3)  is  this 
method  sufficient  for  the  needs  of  SAR;  (4)  is  this  method  acceptable  for 
the  IIP  to  predict  iceberg  motion;  or,  (5)  does  it  need  to  be  corrected  or 
replaced  by  another  model  for  determining  wind  current?  An  experimental 
investigation  is  planned  that  is  designed  to  compare  the  predictions  of  the 
SAR-Mooney  model  with  field  measurements.  The  results  will  be  presented 
in  a  subsequent  report. 


2  Background 

Before  the  merits  of  the  model  can  be  discussed,  it  is  necessary  to  understand 
the  mechanisms  by  which  wind  generates  local  current.  When  wind  blows 
over  the  earth’s  surface,  a  stress  is  exerted  on  the  surface,  whether  this  be 
solid  ground  or  sea.  This  stress  is  a  considerable  retarding  force  on  the 
atmosphere  and  a  driving  force  for  the  ocean.  The  wind  stress  produces  a 
direct  response  called  the  Ekman  component  of  the  current.  This  current 
is  mainly  confined  to  a  thin  layer  (on  the  order  of  10  to  100  m  deep)  near 
the  ocean  surface,  usually  within  the  upper  mixed  layer  (Gill,  1982).  Pure 
drift  currents  are  generated  by  the  drag  of  the  wind  passing  over  homogenous 
water.  At  first,  the  resulting  current  follows  the  direction  of  the  wind,  but  as 
a  consequence  of  the  earth’s  rotation,  the  Coriolis  effect  causes  the  current 
to  deviate  from  the  wind  direction.  As  each  water  layer  moves,  it  exerts  a 
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stress  on  the  layer  beneath  it,  so  that  each  layer  deviates  even  more  from  the 
wind  direction  than  the  layer  above  it. 

Ekman  (1905)  was  the  first  to  examine  the  oceanic  response  to  time- 
dependent  wind  forcing.  By  assuming  an  infinitely  deep  ocean  (to  neglect 
bottom  friction),  unbounded  in  the  horizontal  direction,  with  a  steady  wind 
blowing  over  a  homogenous  ocean  with  a  level  sea  surface  (so  the  pressure 
gradient  is  zero),  the  magnitude  and  direction  of  the  Ekman  velocity  can  be 
found  as  a  balance  between  the  local  accelerations  and  the  Coriolis  force  with 
the  vertical  gradient  of  the  horizontal  stress.  Thus,  the  linearized  equation 
can  be  written: 

dw  d2w 

(i) 

where  w  =  u  +  iv,  the  complex  current,  of  which  u  is  eastward  velocity  and  v 
is  northward  velocity,  /  =  217  sin  <f>  is  the  Coriolis  parameter,  (17  is  the  earth’s 
angular  velocity  and  <f>  is  the  latitude),  v  is  the  kinematic  eddy  viscosity,  and 
z  is  the  vertical  coordinate,  positive  upward.  Only  vertical  derivatives  of 
horizontal  stresses  are  included  in  Equation  1  because  the  vertical  scale  of 
oceanic  boundary  layers  is  much  smaller  than  the  horizontal  scale  over  which 
the  stresses  vary  (Gill,  1982). 

Ekman  assumed  that  the  kinematic  eddy  viscosity  was  constant  with 
depth  to  simplify  the  equations  and  also  because  little  was  known  about 
the  variation  of  the  eddy  viscosity  with  depth.  If  we  assume  that  the  sea 
is  initially  at  rest,  and  that  the  wind  rises  suddenly  and  is  maintained  at  a 
constant  speed,  then  the  solution  to  Equation  1,  assuming  a  constant  eddy 
viscosity,  is  given  by 


d2w 


>  dz 2’ 


(2) 


Jelesnianski  (1970)  obtained  a  solution  i,o  Equation  1  for  the  more  com¬ 
plicated  problem  of  an  arbitrary  time-dependent  wind: 


o  =  jj  £  ^  [(« + jf  -  r)e  Snrdt 


(3) 


where  H  is  the  water  depth,  and 


0n  =  if  + 


(4) 
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with  the  boundary  conditions, 

where  F(t)  is  the  wind  stress.  The  derivation  of  equation  3  can  be  found  in 
Appendix  1. 

Equation  3  forms  the  basis  to  the  manual  method  of  solution  derived  by 
Mooney  (1978). 

3  The  SAR  Manual  Method 

Jelesnianski  (1970)  based  his  model  on  the  assumptions  made  by  Ekman 
(1905),  which  are  listed  below: 

1.  no  horizontal  boundaries  (neglect  boundary  effects); 

2.  infinitely  deep  water  (to  avoid  bottom  friction); 

3.  constant  kinematic  eddy  viscosity  (for  simplification  and  also  be¬ 
cause  little  was  known  about  the  variation  of  the  eddy  viscosity 
with  depth); 

4.  steady  wind  blowing  for  a  long  time  (steady-state  conditions); 

5.  homogeneous  water  and  level  sea  surface  (hydrostatic  pressure); 

6.  constant  Coriolis  parameter. 

Although  these  assumptions  are  too  restrictive  for  problems  with  large  areal 
extent  where  the  variation  in  the  Coriolis  parameter,  density  stratification, 
and  lateral  stresses  are  too  important  to  be  neglected,  Ekman ’s  theory  may 
be  sufficient  for  problems  of  limited  area  where  the  lateral  variation  in  stress 
is  much  smaller  than  the  vertical  variations  (Jelesnianski,  1970).  Verifica¬ 
tion  of  Ekman’s  theory  is  difficult  because  pure  drift  currents  that  satisfy 
the  assumptions  above  are  rarely  found  in  nature  and  because  the  vertical 
distribution  of  the  eddy  viscosity  is  poorly  understood,  thus  precluding  a 
direct  comparison  of  the  theory  with  field  measurements.  However,  it  has 
been  observed  that  simple  models  describe  well  measurements  of  velocity 
fluctuations  in  the  mixed  layer  (Gill,  1982). 
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3.1  Assumptions 

In  addition  to  Ekman’s  simplifying  assumptions,  Mooney  (1978)  made  sev¬ 
eral  other  assumptions  in  deriving  and  tuning  his  model.  The  validity  of  these 
assumptions  and  their  effects  on  the  calculated  wind  drift  current  were  ex¬ 
amined  by  comparison  of  Mooney’s  method  with  the  one- dimensional  version 
of  the  mixed  layer  model  developed  by  Mellor  and  Yamada  (1982).  A  com¬ 
parison  of  the  Mellor- Yamada  Level  2|  model  with  other  commonly  used 
oceanic  mixed-layer  models  (namely,  Mellor- Yamada  Level  2,  Niiler,  Gar¬ 
wood,  Price,  and  Therry-Lacarrere)  can  be  found  in  Martin  (1986).  Com¬ 
parisons  were  made  for  cases  with  both  idealized  forcing  and  data  from  three 
field  experiments.  Noticeable  discrepancies  among  the  models  were  observed. 
In  particular,  significant  differences  were  discovered  for  the  increase  in  the 
depth  of  the  mixed  layer  in  autumn  and  winter  due  to  wind  mixing  and  con¬ 
vection,  and  for  the  decrease  in  depth  of  the  mixed  layer  during  light  winds 
and  strong  heating.  The  Mellor- Yamada  Level  2|  model  was  found  to  un¬ 
derpredict  the  mixed-layer  depth  when  compared  with  field  measurements. 
Additionally,  the  predicted  mixed  layer  was  not  as  well  mixed  as  indicated  by 
the  measurements.  However,  the  model  did  predict  well  the  smooth  gradient 
at  the  base  of  the  mixed  layer  (Martin,  1986). 

3.1.1  Form  of  the  Wind  Stress 

To  solve  Equation  3,  the  form  of  the  wind  stress  must  be  known  or  predicted. 
The  wind  stress  is  usually  expressed  in  terms  of  a  drag  coefficient  and  the 
wind  speed.  The  value  of  the  drag  coefficient  depends  on  the  roughness  of  the 
underlying  surface.  For  winds  over  water,  the  roughness  of  the  water  surface 
can  vary  with  wind  speed  (i.e.,  surface  waves).  Therefore,  the  drag  coefficient 
should  not  be  considered  a  constant,  but  rather  dependent  on  wave  speed. 
Owing  to  the  difficulties  in  wind  stress  estimates,  many  oceanographers  use 
the  quadratic  law,  first  applied  by  Ekman, 

F(t)  =  PaCDW 2  (6) 

where  pa  is  the  density  of  air  and  the  drag  coefficient,  Cd,  is  assumed  equal 
to  2.6xl0-3,  and  W  =  Wx  +  iWy,  is  the  complex  wind  velocity,  Wx  and  Wv 
represent  the  westerly  and  southerly  wind  velocities,  respectively.  Mooney 
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Wind  Speed 
m/s 

Mellor- Yamada 

SAR-Mooney 

d 

km 

Current 

m/s 

Direction 

Current 

m/s 

Direction 

2.5 

0.03 

40 

0.04 

48 

0.35 

5.0 

0.07 

40 

0.08 

48 

0.44 

10.0 

0.16 

30 

0.15 

48 

1.11 

20.0 

0.37 

27 

0.30 

48 

3.11 

40.0 

0.82 

25 

0.60 

48 

8.04 

Table  1 :  Comparison  of  wind  current  and  direction  for  Mellor- Yamada  model 
with  SAR-Mooney  model  at  45°  N.  Water  depth  is  250  m.  Winds  from  the 
south  at  constant  velocity  for  48  hours.  Mellor- Yamada  current  predictions 
at  depth  0.6  m.  The  difference  over  six  hours  in  the  model  predictions,  d ,  is 
given  by  Equation  8. 


assumed  that  the  wind  stress  was  linearly  proportional  to  the  wind  velocity 

F(t)  =  paCDW(t)  (7) 

and  that  the  wind  drag  coefficient,  Co  —  1.9xl0-3  cm/sec. 

Stress  is  expressed  dimensionally  as  mass/length/time2.  So  although  the 
drag  coefficient  is  usually  taken  to  be  dimensionless  Mooney  had  to  assume  a 
dimensional  drag  coefficient  to  make  his  expression  dimensionally  correct.  To 
simplify  the  computations,  Mooney  assumed  a  linear  relation  between  wind 
speed  and  stress,  which  may  restrict  the  method’s  applicability  to  a  limited 
range  of  wind  velocities.  Unfortunately,  Mooney  (1978)  does  not  provide 
any  statistical  information  on  the  wind  data  used  to  tune  his  model,  so  the 
range  of  wind  velocities  over  which  this  method  should  be  most  accurate  is 
not  known.  Therefore,  the  wind  drift  calculated  from  Mooney’s  model  were 
compared  with  the  predictions  from  the  Mellor- Yamada  model  for  a  range 
of  wind  velocities.  Table  1  lists  the  results  for  both  models  for  a  constant, 
southerly  wind  at  45°  N. 

Within  the  range  of  wind  speeds  examined,  the  SAR-Mooney  model  over- 
predicted  the  amount  of  rotation  of  the  current  from  the  wind  direction.  At 
low  wind  velocities  (<  10  m/s  [20  knots]),  this  was  not  significant,  however, 
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the  difference  between  the  current  directions  predicted  by  the  two  models 
increases  with  increasing  wind  velocity.  There  is  a  large  jump  in  the  dis¬ 
agreement  between  the  models  as  the  wind  increases  from  5  to  10  m/s  (10 
to  20  knots),  and  then  a  leveling  off  with  increasing  wind  speed.  The  per¬ 
cent  error  is  greatest  at  the  lowest  wind  speeds.  However,  at  low  currents 
the  distance  drifted  is  less,  so  the  difference  between  the  spatial  drift  over 
time  predicted  by  the  models  is  much  lower  than  at  high  wind  speeds.  The 
difference  in  the  distance  drifted  over  six  hours  can  be  calculated  from  the 
solution  for  oblique  triangles: 

d=  6  [t v2m  +  wlfY  -  2wmwmy  cos  (4>m  —  (!>my)\  1  (8) 

where  the  subscripts  M  and  MY  denote  the  SAR-Mooney  and  Mellor- 
Yamada  predictions,  respectively.  At  10  m/s  (20  knots),  the  SAR-Mooney 
model  changes  from  overpredicting  the  current  to  underpredicting  it  for 
higher  velocity  winds.  However,  before  drawing  conclusions  based  on  these 
results  it  should  be  understood  that  only  one  wind  direction  and  latitude  was 
examined. 

3.1.2  Effect  of  Tuning  Parameters 

The  validity  of  Mooney’s  model  may  be  limited  somewhat  by  the  parameters 
he  used  to  tune  his  model.  Mooney  used  subsurface  current  measurements 
( H  =  — 12  m)  from  approximately  39°  N,  in  water  of  depth  2682  m.  Without 
justification,  Mooney  assumed  that  48  hours  of  wind  speed  and  direction 
data  would  be  sufficient  to  estimate  the  wind  generated  current.  Using  the 
Mellor-Yamada  model  for  comparison,  the  effect  of  tuning  the  model  with 
the  assumed  parameters  are  considered  below.  These  parameters  include  (1) 
the  subsurface  current  record;  (2)  the  water  depth;  (3)  the  length  of  the 
wind  and  current  record;  and  (4)  variations  in  latitude,  in  particular  wind 
drift  predictions  at  the  equator. 

Subsurface  current  record  Mooney  (1978)  tuned  his  model  with  a  time 
series  of  current  measurements  recorded  at  a  depth  of  12  m.  Since  the  effects 
of  the  Ekman  velocity  are  dissipated  with  depth,  it  is  reasonable  to  suspect 
that  the  model  would  underpredict  the  magnitude  and  rotation  of  the  sur¬ 
face  current.  However,  when  compared  with  the  Mellor-Yamada  model,  the 
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Wind  Speed 
m/s 

Mellor-Yamada 

SAR-Mooney 

d 

km 

Current 

m/s 

Direction 

Current 

m/s 

Direction 

2.5 

0.01 

76 

0.04 

48 

0.70 

5.0 

0.03 

76 

0.08 

48 

1.20 

10.0 

0.08 

59 

0.15 

48 

1.72 

20.0 

0.20 

49 

0.30 

48 

2.33 

40.0 

0.47 

43 

0.60 

48 

3.07 

Table  2:  Comparison  of  wind  current  and  direction  for  Mellor-Yamada  model 
with  SAR-Mooney  model  at  45°  N.  Water  depth  is  250  m.  Winds  from  the 
south  at  constant  velocity  for  48  hours.  Mellor-  Yamada  current  predictions 
at  depth  6.9  m.  The  difference  over  six  hours  in  the  model  predictions,  d,  is 
given  by  Equation  8. 

Mooney  model  predictions  were  inconsistent.  Table  2  lists  the  currents  at 
a  depth  of  6.9  m  predicted  by  the  Mellor-Yamada  model.  For  wind  speeds 
under  20  m/s  (40  knots),  Mooney’s  model  underpredicted  the  rotation  of  the 
current  in  comparison  with  Mellor’s  model  and  for  the  range  of  wind  veloci¬ 
ties  examined,  Mooney’s  model  overpredicted  the  magnitude  of  the  current. 
Mooney’s  predictions  were  in  better  agreement  with  Mellor’s  near  the  surface 
(see  Table  1)  than  near  the  depth  of  his  current  data  which  he  used  to  de¬ 
termine  the  coefficients.  For  wind  drift  calculations,  we  are  not  particularly 
interested  in  the  currents  except  in  the  vicinity  of  the  surface,  however,  it 
should  be  noted  that  the  observed  better  agreement  near  the  surface  of  the 
predictions  is  serendipitous. 

Water  depth  The  SAR  manual  states  that  wind  current  should  be  deter¬ 
mined  for  water  depths  greater  than  30  m  (100  ft).  Since  the  coefficients 
for  Mooney’s  model  were  calculated  in  much  deeper  water  and  one  of  the 
assumptions  in  Ekman’s  theory  is  that  the  water  is  infinitely  deep,  it  is  rele¬ 
vant  to  ascertain  if  this  model  is  valid  in  shallower  water.  Ekman  assumed  an 
infinitely  deep  ocean  so  the  surface  and  the  bottom  frictional  layers  (known 
as  the  Ekman  layers)  would  not  meet.  Where  the  water  depth  is  on  the 
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order  of  the  Ekman  depth  or  less,  the  surface  and  bottom  layer  meet  and 
can  even  overlap.  The  two  spirals  tend  to  cancel  each  other  with  the  result 
that  the  total  transport  is  more  in  the  direction  of  the  surface  wind.  When 
the  water  depth  decreases  to  about  one  tenth  of  the  Ekman  depth,  the  trans¬ 
port  is  essentially  in  the  direction  of  the  wind  because  the  Coriolis  effect  is 
overpowered  by  frictional  effects  (Pond  and  Pickard,  1983). 

The  Ekman  depth, 

D  =  (9) 

where  p  is  the  density  of  sea  water,  is  independent  of  the  wind  stress,  however, 
it  does  increase  with  increasing  eddy  viscosity  and  with  decreasing  latitude. 
At  the  equator,  the  Ekman  depth  approaches  infinity.  Therefore,  the  water 
depth  will  become  increasingly  more  important  as  the  equator  is  approached. 

As  the  wind  speed  increases,  the  eddy  viscosity  increases,  so  the  Ekman 
depth  increases.  Therefore  at  high  wind  speeds  and  approaching  the  equator, 
30  m  is  insufficient  for  the  Ekman  spiral  to  develop.  An  absolute  limit  for  the 
applicability  of  the  model  is  not  practical  due  to  the  complex  relationship 
between  wind  speed,  latitude,  water  depth  and  the  Ekman  depth.  Figure  1 
shows  a  plot  that  can  be  used  to  estimate  the  water  depth  limit  for  which  the 
model  is  valid.  Use  of  the  SAE  model  in  water  depths  shallower  than  shown 
in  Figure  1  may  lead  to  incorrect  estimates  of  drift  velocity  and  direction. 

When  the  ratio,  H/D,  between  the  water  depth,  H,  and  the  Ekman 
depth,  D,  is  small,  the  angle  between  the  wind  and  the  current  direction, 
known  as  the  deflection  angle,  7,  is  small,  and  the  surface  current  flows  nearly 
in  the  direction  of  the  wind.  As  H/D  increases,  7  is  alternately  smaller  and 
larger  than  45°.  When  the  oceanic  depth  is  approximately  equal  to  the  Ek¬ 
man  depth,  the  assumption  of  infinitely  deep  water  is  reasonable  (Neumann 
and  Pierson,  1966).  The  following  table  from  Neumann  and  Pierson  shows 
deflection  angles  for  different  depth  ratios: 


Thtd 

0.25 

0.5 

0.75 

0 

0 

H 

(KH 

21.5 

45 

45.5 

45 

Length  of  wind  and  current  record  The  SAR-Mooney  method  uses 
a  wind  data  interval  of  6  hours  for  practical  rather  than  physical  reasons; 
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Figure  1:  Ekman  depth  as  a  function  of  wind  speed  and  latitude. 
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the  synoptic  weather  charts  are  for  six  hour  periods.  Without  explanation, 
Mooney  used  a  time  record  length  of  48  hours,  which,  divided  by  6-hour 
intervals,  results  in  eight  time  periods.  The  use  of  a  48  hour  record  length 
without  justification  is  a  matter  of  concern. 

Ekman’s  theory  assumes  an  ocean  initially  at  rest;  however,  this  is  rarely 
found  in  nature.  Inertial  oscillations  in  the  ocean  often  grow  to  large  am¬ 
plitude  and  then  decay  slowly,  therefore,  the  length  of  the  wind  record  must 
be  long  enough  for  wind  generated  currents  developed  by  pre-existing  wind 
conditions  to  have  sufficient  time  to  decay.  Studies  indicate  a  decay  time  of 
several  inertial  periods  (Gill,  1982)  where  the  inertial  period  is 

Ti  =  2*7/.  (10) 


Figures  2a  and  b  show  the  inertial  period  as  a  function  of  latitude.  Since  iner¬ 
tial  oscillations  decay  more  slowly  as  the  equator  is  approached,  the  possibil¬ 
ity  of  the  current  estimates  being  incorrect  due  to  the  presence  of  pre-existing 
inertial  oscillations  increases. 

A  sudden  change  in  wind  can  also  reduce  or  even  eliminate  pre-existing 
oscillations.  Assume  inertial  oscillations  are  present  when  the  wind  changes 
at  time  t  =  t0,  and  w0  is  the  value  of  the  Ekman  current  at  time  t  —  t0,  then 
the  total  Ekman  current  is, 


w 


—ii>  d2w  ( 
f  dz 2  +  1^°  + 


iv  d2w\ 
7~dz2  / 


g  —  •’/(*- *o) 


(11) 


where  the  amplitude  of  the  inertial  oscillations  is  given  by  the  factor  in  curly 
brackets.  If  the  wind  change  occurs  at  a  phase  of  the  pre-existing  oscillations 
when  this  factor  is  small,  the  resulting  oscillations  will  be  weak.  If  the 
wind  change  occurs  when  this  factor  is  large,  strong  oscillations  will  occur 
(Gill,  1982).  The  effect  of  the  time  of  the  wind  change  on  the  amplitude  of 
the  inertial  oscillations  is  shown  in  Figure  3. 

It  may  not  be  feasible  in  all  cases  (particularly  in  the  equatorial  region)  to 
go  back  sufficiently  far  in  the  wind  record  to  allow  for  the  decay  of  any  pre¬ 
existing  inertial  oscillations.  But  for  cases  where  a  strong  storm  has  passed 
recently  through  the  region,  the  possibility  of  pre-existing  currents  should  be 
considered. 


Variations  in  latitude  Table  3  compares  the  predictions  of  the  Mellor- 
Yamada  model  with  those  of  the  SAR-Mooney  model  for  several  latitudes. 
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Figure  2:  Variation  of  inertial  period  with  latitude;  (a)  latitude  1°  to  5°,  (b) 
latitude  5°  to  80°. 
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Figure  3:  A  progressive  vector  diagram  for  a  particle  when  a  westerly  wind 
suddenly  rises  at  time  t  =  0.  Later,  the  wind  changes  to  the  northwest  and 
the  wind  stress  doubles.  The  solid  line  depicts  the  result  when  the  change 
occurs  at  5/4  inertial  periods  after  t  =  0.  The  dashed  line  shows  the  result 
if  instead,  the  wind  change  occurs  at  3/4  of  am  inertial  period  after  t  =  0 
(redrawn  from  Gill,  1982).  For  further  discussion,  see  Gill,  (1982). 
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Latitude 
o  N 

Mellor-Yamada 

SAR-Mooney 

d 

km 

Current 

m/s 

Direction 

Current 

m/s 

Direction 

5 

0.35 

29 

0.68 

40 

7.30 

45 

0.08 

59 

0.15 

48 

1.72 

65 

0.16 

35 

0.15 

54 

1.13 

Table  3:  Comparison  of  wind  current  and  direction  for  Mellor-Yamada  model 
with  SAR-Mooney  model  for  a  range  of  latitudes.  Water  depth  is  250  m. 
Winds  from  the  south  at  10  m/s  (20  knots)  for  48  hours.  Mellor-Yamada 
current  predictions  at  depth  —6.9  m. 

The  difference  between  the  two  models  is  much  greater  in  the  vicinity  of 
the  equator.  At  65°  N,  the  agreement  is  close  for  the  current  speed,  but 
less  good  for  the  current  direction.  Further  south  at  45°  N,  the  agreement 
of  the  current  speeds  is  less  good  than  at  65°  N,  but  the  directions  are  in 
better  agreement,  and  near  the  equator,  at  5°  N,  the  agreement  between 
the  velocities  is  even  worse,  though  the  difference  between  the  predicted 
directions  has  only  increased  slightly.  This  seems  to  indicate  that  the  errors 
incurred  using  the  SAR-Mooney  model  increase  as  the  equator  is  approached. 

3.2  Additional  Considerations 

3.2.1  Drift  near  the  equator 

In  the  vicinity  of  the  equator,  the  inertial  frequency  goes  to  zero  (/  = 
2f l  sin  <f>  =  0)  so  that  the  wind  drift  is  dependent  solely  on  the  wind  stress, 
and  therefore,  the  drift  will  essentially  be  in  the  direction  of  the  wind. 

3.2.2  Validity  of  model  in  nearshore  region 

The  SAR  manual  requires  the  wind  current  to  be  calculated  when  the  dis¬ 
tance  from  shore  is  greater  than  32  km  (20  miles).  This  is  based  on  the 
effect  of  a  coastal  boundary  on  the  development  of  wind  generated  currents. 
Ekman  assumed  an  infinite  ocean  so  that  the  effects  of  boundaries  could  be 
neglected. 
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In  the  vicinity  of  a  barrier,  such  as  a  coastline,  a  convergence  or  divergence 
of  the  Ekman  flow  occurs,  resulting  in  a  piling  up  of  water  or  in  a  reduction 
of  the  sea  level,  depending  on  how  the  wind  is  directed  with  respect  to  the 
coastline.  In  shallow  water,  flow  toward  the  coast  causes  water  to  pile  up 
at  the  shoreline.  Conversely,  an  offshore  flow  will  cause  the  sea  surface  to 
be  lowered  near  the  coast.  Therefore,  the  wind  direction  with  respect  to  the 
coastline  is  important,  except  where  the  water  depth  is  considerably  smaller 
than  twice  the  Ekman  depth.  In  this  case,  the  Coriolis  effect  is  negligible 
(Neumann  and  Pierson,  1966). 

The  presence  of  the  shoreline  will  effect  the  magnitude  and  direction  of 
a  wind  generated  current  where  the  distance  to  the  coast  is  of  order  of  the 
Rossby  radius: 

«  =  W'Vl/l-  (12) 

Thus,  in  30  m  of  water  (the  SAR  Manual  recommended  depth  limit  for 
using  the  SAR-Mooney  model)  at  45°  N,  the  Rossby  radius  is  about  160  km 
(100  miles).  Within  a  Rossby  radius  of  the  shoreline,  a  longshore  wind  will 
generate  a  current  towards/away  from  the  coast.  Since  there  can  be  no  flow 
across  the  coastline,  this  steady  current  will  cause  a  rise/fall  in  the  sea  level 
at  the  coast  in  order  to  conserve  mass.  Therefore,  the  sea  level  and  the 
longshore  current  will  change  continually  with  time  (Gill,  1982). 

Additionally,  wind  drift  will  only  be  affected  by  the  presence  of  the  coast¬ 
line  when  the  wind  duration  is  sufficiently  long  for  a  surface  gravity  wave  to 
travel  from  the  shoreline  to  the  point  of  interest.  The  time  necessary  can  be 
calculated  from 

(  =  £/(sff),/2  (13) 

where  L  is  the  distance  to  the  shore. 

Changes  in  the  wind  field  can  affect  the  Ekman  transport  in  a  similar 
manner  as  a  coastal  boundary.  For  instance,  consider  the  North  Pacific 
or  North  Atlantic  around  25°  N  between  the  prevailing  westerlies  and  the 
Northern  Trade  Winds.  At  the  boundary  between  the  wind  directions,  the 
Ekman  transport  convergences  because  both  winds  transport  the  water  to 
the  right  of  the  wind.  This  also  occurs  in  the  Southern  Hemisphere  with  a 
transport  to  the  left  (Neumann  and  Pierson,  1966). 
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3.3  Summary 

The  National  SAR  Manual  method  for  calculating  wind  driven  currents  has 
been  studied  to  determine  its  usefulness  in  search  and  rescue  applications. 
The  derivations  of  the  model  have  been  reviewed  to  verify  that  it  is  physically 
correct.  It  has  been  compared  with  one  of  the  most  complex  one  dimensional 
models  of  the  mixed  layer  currently  available  in  the  literature.  Several  aspects 
of  the  model  have  been  assessed  and  the  conclusions  are  summarized  below. 

3.3.1  Form  of  the  wind  stress 

In  comparison  with  the  Mellor-Yamada  model,  the  SAR-Mooney  model  over¬ 
estimated  the  deviation  of  the  current  direction  from  the  wind  direction  for 
the  entire  range  of  wind  velocities  examined  (2.5  to  40  m/s  [5  to  80  knots]). 
The  model  predictions  of  current  velocity  agreed  well  with  a  steady  wind 
speed  of  10  m/s  (20  knots),  but  the  difference  between  the  models  increased 
with  increasing  windspeed. 

3.3.2  Tuning  parameters 

Subsurface  current  record  The  SAR-Mooney  model  was  tuned  with 
subsurface  current  data  from  12  m,  however,  when  compared  with  the  Mellor- 
Yamada  model,  it  agreed  with  near-surface  predictions.  Although  this  result 
suggests  that  the  SAR-Mooney  model  may  not  describe  well  the  experimental 
data  on  which  its  coefficients  are  based,  it  does  not  appear  to  be  a  drawback 
to  its  use. 

Water  depth  Ekman’s  theory  assumes  infinitely  deep  water  and  Mooney 
(1978)  used  data  from  relatively  deep  water  compared  with  the  depth  limita¬ 
tion  recommended  by  the  National  SAR  Manual  of  30  m  (100  feet).  Where 
the  water  depth  is  equal  to  or  greater  than  the  Ekman  depth,  given  by  Equa¬ 
tion  9,  the  wind  drift  should  not  be  affected  by  the  water  depth.  At  depths 
shallower  than  the  Ekman  depth,  the  SAR-Mooney  predictions  should  be 
treated  with  caution.  Figure  1  may  be  used  to  estimate  the  limiting  depth 
for  which  the  water  can  be  assumed  to  be  infinitely  deep  as  a  function  of 
windspeed  and  latitude. 
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Variations  in  latitude  Mooney  (1978)  tuned  his  model  with  wind  and 
current  data  from  about  39°  N.  It  appears  that  the  SAR-Mooney  model  may 
be  increasingly  inaccurate  as  the  equator  is  approached,  particularly  with 
respect  to  the  magnitude  of  the  current. 


Equatorial  region  In  the  vicinity  of  the  equator,  the  inertial  frequency 
approaches  zero,  so  the  wind  drift  is  dependent  only  on  the  wind  stress  and 
Equation  1  reduces  to 


dw  __ 

~dt  =  V!h* 

so  the  drift  will  be  in  the  direction  of  the  wind. 


(14) 


Nearshore  region  The  recommended  limitation  of  the  SAR-Mooney 
method  to  cases  greater  than  32  km  (20  miles)  from  the  coast  appears  to 
be  too  short  a  distance.  Wind  drift  will  be  affected  by  a  coastal  barrier  at 
a  distance  equal  to  or  less  than  the  Rossby  radius,  (Equation  12).  As  the 
water  depth  decreases,  the  effect  of  the  coastline  increases.  Additionally,  the 
effect  of  the  coast  increases  linearly  with  decreasing  distance  from  the  coast. 


4  Verification  by  Field  Measurements 

The  investigation  described  in  this  report  was  limited  to  the  comparison  of 
two  models  of  wind  generated  current;  one  was  a  complicated  one  dimen¬ 
sional  mixed-layer  model  while  the  other  was  a  simple  model  dependent  on 
carefully  determined  coefficients.  It  is  useful  to  compare  these  models  not 
just  against  each  other,  but  also  with  experimental  data.  An  experiment 
is  planned,  designed  to  test  the  accuracy  of  the  SAR-Mooney  model.  The 
measurements  obtained  from  the  experiment  will  be  presented  in  a  future 
report  with  the  results  of  the  model.  The  comparison  of  the  measurements 
and  model  predictions  will  be  used  to  ascertain  the  usefulness  and  limitations 
of  the  model  for  SAR  applications. 
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5 


Recommendations 


The  following  recommendations  are  based  on  the  comparison  of  the  predic¬ 
tions  of  the  SAR-Mooney  model  with  those  of  the  Mellor-Yamada  Level  2| 
model.  Consequently,  these  recommendations  should  be  regarded  as  provi¬ 
sional  until  the  experimental  investigation  has  been  performed  and  the  results 
analyzed. 

5.1  Manual  Method 

During  this  interim  period,  the  SAR-Mooney  model  should  continue  to  be 
used  for  search  and  rescue  and  iceberg  drift  applications  with  the  following 
caveats: 

1.  The  SAR  model  over-estimates  the  deviation  of  the  current  from 
the  wind  direction. 

2.  At  high  wind  speeds  (over  10  m/s  [20  knots]),  the  SAR  model 
over-estimates  the  current  speed. 

3.  At  depths  shallower  than  the  Ekman  depth,  the  results  should 
be  treated  with  caution.  Figure  1  may  be  used  to  estimate  the 
limiting  depth  for  which  the  water  can  be  assumed  to  be  infinitely 
deep  as  a  function  of  windspeed  and  latitude. 

4.  Although  at  mid  and  high  latitudes,  the  model  predicts  well  the 
current  speed,  the  SAR-Mooney  model  may  be  increasingly  inac¬ 
curate  as  the  equator  is  approached,  particularly  with  respect  to 
the  magnitude  of  the  current. 

5.  In  the  vicinity  of  the  equator,  the  drift  will  be  in  the  direction  of 
the  wind. 

6.  The  recommended  limitation  of  the  SAR-Mooney  method  to  cases 
greater  than  32  km  (20  miles)  from  the  coast  appears  to  be  too 
short  a  distance.  Wind  drift  will  be  affected  by  a  coastal  barrier  at 
a  distance  equal  to  or  less  than  the  Rossby  radius,  a  =  (gH)1^2 /\f\. 

7.  As  the  water  depth  decreases,  the  effect  of  the  coastline  increases. 
Additionally,  the  effect  of  the  coast  increases  linearly  with  decreas¬ 
ing  distance  from  the  coast. 
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5.2  Numerical  Model 


The  SAR- Mooney  model  should  be  replaced  in  CASP  by  a  more  complex  one 
dimensional  mixed  layer  model.  Several  models  are  discussed  in  the  litera¬ 
ture.  These  models  differ  in  terms  of  their  formulation,  complexity,  and  effi¬ 
ciency  (Martin,  1986).  The  models  can  generally  be  described  as  differential 
or  bulk  models,  depending  on  their  method  of  formulation.  The  differential 
models  use  the  equations  for  the  conservation  of  momentum,  heat,  salt,  and 
turbulent  kinetic  energy  in  their  primitive  form  and  are  not  integrated  over 
the  mixed  layer.  The  bulk  models  assume  that  the  mixed  layer  is  well-defined 
and  the  equations  of  conservation  are  integrated  over  the  mixed  layer  (Mar¬ 
tin,  1986).  The  most  appropriate  model  for  incorporation  into  CASP  will 
depend  on  the  model  formulation  and  computation  ease. 
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A  Appendix 


The  linearized  equation  of  motion  for  a  pure  drift  current  is  given  by  Ek- 
man  (1905)  as: 

dw  d2w 

—  =  (15) 

where  w  =  u  +  iv  is  the  complex  current  of  which  u  is  the  eastward  velocity 
and  v  is  the  northward  velocity.  The  Coriolis  parameter  is  given  by,  /  = 
2ft  sin  <f>  where  ft  is  the  earth’s  rotation  rate  and  phi  is  the  latitude.  The 
kinematic  eddy  viscosity  is  represented  by  v  and  z  is  the  vertical  coordinate, 
positive  upward. 

The  initial  conditions  are  at  rest: 

r\ 

w(z,0)  =  -£■  =0,  (16) 

91  (=0 

and  boundary  conditions  are  no  slip  at  the  ocean  bed: 

=  0  (17) 

where  the  water  depth  is  given  by  H,  and  the  wind  stress  at  the  ocean  surface, 
F,  is  a  function  of  time: 

*  ^  =  F(t).  (18) 

uz  1=0 

The  method  of  solution  to  Equation  15  for  an  arbitrary  and  unsteady 
wind  follows  that  given  by  Jelesnianski  (1970).  Equation  15  can  be  solved 
using  a  Laplace  transform,  such  that: 


su>(z,s)  =  — ifw(z,s )  +  v 


with  boundary  conditions 


w(—H,s)  =  0,  v 


d*w(z,  s) 


=  Hs). 


Again,  following  Jelesnianski  (1970),  assume  the  surface  stress  is  am  instan¬ 
taneously  applied  constant  force  iF0  with  the  resulting  current,  w0.  Then, 
Equation  19  can  be  solved 

iF0  sinha(/f  +  z)  .  x 

w°  = - TIT -  (21) 

vo.  cosh  Ha 


A-l 


where 


a  =  ^~ 


if  +  s 


(22) 


Equation  21  is  solved  using  the  inverse  Laplacian  and  applying  the  s-shift 
property; 

+  a))  =  e-"L-'(f(*))  =  e~atf(t).  (23) 

The  inverse  Laplacian  can  then  be  found  as 


w0  =  ^E~o(-l)nsin[(n+|)(/f  +  ^] 
exp  +  +  *]  . 


(24) 


(see  for  instance,  Oberhettinger  and  Badii,  1973). 

If  the  restriction  of  a  constant  wind  stress  is  replaced  with  a  temporally 
varying  wind,  initially  at  rest,  the  solution  to  Equation  19  can  then  be  given 

by 

F(s)sinha(/f  +  *)  F(s)w0 

w{z,s)  = - r- 77 - =  — r= — .  (25) 

ua  cosh  Ha  iF0 

Solving,  using  Laplace  convolution, 


1  f* 

w  =  — -  /  Fit  -  T)w0dr. 
iF0  Jo 


(26) 


Inserting  Equation  24  into  Equation  26,  the  wind  drift  current  can  be  found 
from 

w(z ,  t)  =  -jj  £  COS  [(n  +  Jq  F(t  -  T)e~KTdr  (27) 

where  0n  is  given  by 

#„=;/+[  <28> 


A-2 


